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ABSTRACT 

Cells release RNA-carrying vesicles and membrane- 
free RNA/protein complexes into the extracellular 
milieu. Horizontal vesicle-mediated transfer of 
such shuttle RNA between cells allows dissemin- 
ation of genetically encoded messages, which may 
modify the function of target cells. Other studies 
used array analysis to establish the presence of 
microRNAs and mRNA in cell-derived vesicles 
from many sources. Here, we used an unbiased 
approach by deep sequencing of small RNA 
released by immune cells. We found a large variety 
of small non-coding RNA species representing per- 
vasive transcripts or RNA cleavage products 
overlapping with protein coding regions, repeat se- 
quences or structural RNAs. Many of these RNAs 
were enriched relative to cellular RNA, indicating 
that cells destine specific RNAs for extracellular 
release. Among the most abundant small RNAs in 
shuttle RNA were sequences derived from vault 
RNA, Y-RNA and specific tRNAs. Many of the 
highly abundant small non-coding transcripts in 
shuttle RNA are evolutionary well-conserved and 
have previously been associated to gene regulatory 
functions. These findings allude to a wider range of 
biological effects that could be mediated by shuttle 
RNA than previously expected. Moreover, the data 
present leads for unraveling how cells modify the 
function of other cells via transfer of specific 
non-coding RNA species. 



INTRODUCTION 

Nano-sized membrane vesicles represent a recently 
identified class of intercellular communication vehicles 
operating in many organisms (1-6). Such vesicles can 
derive from multivesicular bodies (MVBs), which are 
late endosomal compartments containing multiple 
50-100 nm sized intraluminal vesicles. Fusion of MVBs 
with the plasma membrane causes the release of their 
intraluminal vesicles, which are then called exosomes (7). 
Alternatively, vesicles can be released by cells through 
direct shedding from the plasma membrane (1). Cells 
can tightly regulate the release and molecular composition 
of these vesicles (8,9) and vesicle targeting depends on the 
type and activation status of recipient cells (10,11). 
Despite their early description decades ago (12,13), the 
wide-spread occurrence of cell-derived vesicles and their 
potential for tailor-made modulation of target cell 
behaviour has only been recognized during the last few 
years. It is now clear that cell-derived vesicles are not 
only released by almost all cultured cell types, but are 
also present in a wide range of body fluids (1). Since the 
molecular make-up and release of cell-derived vesicles is 
regulated by the producing cell, these vesicles are of 
interest for disease-related biomarkers (14). Moreover, 
extracellular vesicles may be used as therapeutic agents 
(15,16). 

Besides specific sets of lipids and proteins, cells can 
shuttle RNA into vesicles determined for release into the 
extracellular space. This allows the conveyance of genet- 
ically encoded messages between cells (17). The first key 
publication on nucleic acids in cell-derived vesicles 
reported the presence of miRNA and mRNA in vesicles 
derived from mast cells and the functional transfer of 
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RNA to vesicle-targeted cells (17). More recently, the 
luminal protein and RNA contents of cell-derived 
vesicles were demonstrated to be delivered into the cyto- 
plasm of recipient cells via fusion of vesicles with these 
cells (18). It is currently not known whether all vesicle 
populations released by cells contain RNA. Although 
several studies indicated that the RNA composition of 
cell-derived vesicles is different from the parental cell 
(17,19-21), it is unknown how RNAs are selected for se- 
cretion into the extracellular space. Furthermore, extracel- 
lular RNA may also be associated with macromolecular 
complexes that are not enclosed by a vesicle. We here 
collectively refer to extracellular RNA as 'shuttle RNA'. 
Although the majority of circulating miRNAs in human 
plasma and serum were found to co-fractionate with 
proteins such as Argonaute2 (Ago2) and nucleophosmin, 
which have been suggested to protect miRNAs from deg- 
radation in RNase-rich environments (22,23), it is cur- 
rently unknown if vesicle-free RNAs can bind and 
modify the function of target cells. 

Next-Generation Sequencing (NGS) techniques have 
led to the discovery of large numbers of unexpected 
non-coding RNAs (24). These transcripts were found to 
overlap with exons, introns and intergenic regions (24-30). 
Various non-coding transcripts found in close vicinity to 
protein coding genes, such as promoter-associated long 
and short RNAs, transcription start site-associated 
RNAs, and PROMoter uPstream Transcripts 
(PROMPTS) are suspected to act as regulatory elements 
to modulate gene activity (30). Interestingly, many small 
non-coding RNA species have been found that could act 
as regulatory RNAs similar to miRNAs. Fragments 
derived from small nucleolar RNA (snoRNA), vault 
RNA (vRNA) and transfer RNA (tRNA), for example, 
were shown to bind Argonaute (AGO) proteins and form 
RNA-induced silencing complexes (RlSCs) capable of 
regulating expression of target niRNAs analogous to 
miRNA-containing RISCs (31-33). 

The analysis of shuttle RNA populations has almost 
exclusively focused on miRNAs and mRNAs, most 
likely due to the availability and ease of array hybridiza- 
tion techniques to detect these RNAs. Which other small 
RNA biotypes were released by cells via vesicles or in 
protein complexes was unknown. We therefore aimed to 
comprehensively analyze small RNA species released by 
cells using NGS. Among the most weU studied cell-derived 
vesicles are those produced by dendritic ceUs (DCs). These 
cells are the sentinels of the immune system and the inter- 
action of DC with T cells is the key event in initiation of 
immune responses. We previously found that DC-T cell 
interactions induce vesicle release by DC and that 
DC-derived vesicles can be targeted to activated T cells 
(8,11). The release and targeting of RNA-containing 
vesicles by DC and T ceUs may play an important role 
in the communication between these ceUs and the 
ensuing immune response. Here, we extensively 
characterized smaU (<70nt) shuttle RNA released 
during DC-T cell interactions and compared this RNA 
with smaU cellular RNA in order to investigate selectivity 
in the extracellular release of RNA species. 



MATERIALS AND METHODS 

Cell culture 

Dl ceUs were maintained in Iscove's Modified Dulbecco's 
medium supplemented with 2mM Ultraglutamine 
(Biowhittaker), 10% heat inactivated fetal calf serum 
(FCS, Sigma- Aldrich), lOOIU/ml penicillin and lOOmg/ 
ml streptomycin (GIBCO), 50 nM p-mercaptoethanol 
and 30% conditioned medium from GM-CSF producing 
NIH 3T3 ceUs (Rl). The p53-specific CD4+ T-cell clone, 
generated in a C57BL/6 p53-/- mouse (34) and provided 
by Prof C. Melief (Leiden University Medical Center, 
Leiden, The Netherlands), was cultured as described pre- 
viously (8). For cognate DC/T cell co-culture conditions, 
DC were pre-loaded with 2.5 |iM p53 peptide (amino acid 
77-96) for 2 h and subsequently mixed in a 1:1 ratio with T 
cells, and co-cultured for 20 h in medium containing over- 
night ultracentrifuged (100000^) FCS and Rl conditioned 
medium (to deplete bovine and Rl derived vesicles and 
large protein aggregates). All cultures were maintained 
at 37°C, 5% CO2. CeU death in DC and T ceU cultures 
was less than 5% and cognate DC-T ceU interactions did 
not lead to additional cell death in the indicated co-culture 
period. Experiments were approved by the institutional 
ethical animal committee at Utrecht University (Utrecht, 
The Netherlands). 

Shuttle RNA isolation 

15 X 10*' p53 peptide-pulsed DCs were co-cultured with 
15x 10*" T cells in 25 ml culture medium as described 
earlier. Supernatants were pooled to a total of 450 ml. 
This culture supernatant was subjected to differential cen- 
trifugation (35). In short, cells were removed by two se- 
quential centrifugations at 200g for 10 min. Collected 
supernatant was subsequently centrifuged two times at 
500g for 10 min, foUowed by lOOOOg for 30 min. Vesicles 
and protein/RNA complexes were finally pelleted by ultra- 
centrifugation at lOOOOOg for 65 min in SW28 and SW40 
rotors (Beckman). When indicated, the lOOOOOg sedi- 
mented material was further fractionated using gradient 
density centrifugation. To this end, the lOOOOOg sedi- 
mented material was mixed with 2.5 M sucrose, overlaid 
with a hnear sucrose gradient (2.0-0.4 M sucrose in PBS) 
and floated into the gradient by centrifugation in a SW40 
tube (Beckman) for 16 h at 192 OOOg. The bottom two frac- 
tions (1.26-1.28 g/ml, as measured by refractometry) and 
fractions with densities of 1.12-1.18 g/ml were pooled, 
diluted with PBS/0.1% BSA and centrifuged again at 
lOOOOOg for 65 min. SmaU RNA (<200nt) was isolated 
from lOOOOOg pellets and from co-cultured DC and T 
cells (15 X 106 each) using the mirVana^"^ miRNA 
Isolation Kit (Ambion) according to the manufacturer's 
procedure. The RNA integrity and size distributions were 
assessed using Agilent 2100 Bioanalyzer pico-RNA chips. 

Small RNA sequencing library generation 

Sequencing libraries were generated as described (36). In 
short, small RNA fractions were ligated to adaptors using 
the SOLID Small RNA Expression Kit (Ambion). 
Libraries were pre-amplified using primers containing 
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sequences that make the SREK hbraries compatible with 
the lllumina flow ceU. DNA was denatured for 30" at 98°C 
foUowed by 18 cycles with 30" at 98°C, 30" at 65°C, 30" at 
72°C and final extension for 5' at 72°C. Library fragments 
were separated on a native 6% gradient pre-cast PAGE 
gel (Novex, Invitrogen). The 100-150 bp size fractions 
containing inserts of 20-70 nt in length were excised, 
DNA was eluted from the gel, precipitated and dissolved 
in nuclease free water. DNA yield was quantified using an 
Agilent Bioanalyzer high sensitivity DNA chip. Small 
RNA inserts were single end sequenced for 35 cycles 
using standard lllumina protocols. 

Small RNA sequence data analysis 

Adaptor trimming, genome ahgnment and smaU RNA an- 
notation were performed using the E-miR data analysis 
pipehne as described (36). Sequence reads were aligned to 
the Mouse Mm9 reference genome with bowtie, allowing 
for ahgnments to three loci with two mismatches, foUowed 
by a crossmapping correction as described by de Hoon 
et al. (37). For further annotation with coding and 
repeat RNAs, reads with a minimum overlap of 1 nt 
were clustered into regions and regions containing a 
total of more than five reads in shuttle and cellular 
RNA were annotated using custom perl scripts (38) with 
exon, 3'-UTR, 5'-UTR, intron and RNA biotype informa- 
tion from Ensembl biomart (39), all in a strand-specific 
manner. Intersections with repeat regions were determined 
using information from the repeat track in the UCSC 
genome browser (NCBI137/mm9) and the Galaxy inter- 
face, considering a minimum of 1 nt overlap between the 
sequenced and the repeat region. To determine reproduci- 
bility of the sample preparation and sequencing protocols, 
we have evaluated the Pearson's correlation between the 
abundance of small RNAs in replicate experiments. A 
square root transformation on the number of counts per 
smaU RNA was apphed to stabilize variance and reduce 
influence of high-abundant smaU RNAs on the Pearson's 
correlation coefficient (38). 

Quantitative real-time PCR 

Sequencing derived expression profiles were compared 
with qPCR measurements on small RNA fractions from 
independent biological samples. cDNA was generated 
using the Ncode miRNA First-Strand cDNA Synthesis 
kit (Invitrogen). An equivalent of 100 pg RNA was used 
in the PCR reactions with 100 nM primers (Isogen) and 
4 ul SYBR Green Sensimix (Biohne) in an 8 ^l reaction. 
PCR amplification efficiencies were determined using 
10-fold dilution series of template DNA and were 
between 1.8 and 2.3 for all primer sets. Transcript 
specific forward primers for miR-29a, miR-155 and 
miR-191 were based on full-length miRbase sequences 
(40): 

3p mmu-mir-29a forward: 5'-TAGCACCATCTGAAAT 
CGGTTA 

5p mmu-mir-155 forward: 5'-GGGTTAATGCTAATTG 
TGATAGGG 



5p mmu-mir-191 forward: 5'-CGGAATCCCAAAAGCA 
GCTG 

All forward miRNA primers were used in combination 
with the NCode universal reverse primer. 

Y-RNA forward: 5'-GTGTTTACAACTAATTGATCA 
CAACC 

Y-RNA reverse: NCode universal reverse primer 
SRP-RNA forward: 5'-GGAGTTCTGGGCTGTAGT 
GC 

SRP-RNA reverse: 5'-ATCAGCACGGGAGTTTTG 
AC 

Cycling conditions were as follows: 95°C for lOmin 
followed by 50 cycles of 95°C for 10 s, 57°C for 30 s and 
72°C for 20 s. All PCR reactions were performed on the 
Bio-Rad iQ5 Multicolor Real-Time PCR Detection 
System (Bio-Rad). Raw threshold cycle (Ct) values were 
calculated using the iQ5 Optical System software v. 2.0 
using automatic basehne settings. Thresholds were set in 
the hnear phase of the amplification plots. 



RESULTS 

Sample preparation and sequencing 

Specific antigen-bearing DCs were co-cultured with 
cognate T cells, leading to mutual activation of both cell 
types. Vesicles and large protein complexes released 
during these interactions were isolated from cell culture 
supernatant by differential centrifugation, with a final pel- 
leting step at 100000^ (35). Total RNA was isolated from 
this fraction and we observed that the majority of this 
shuttle RNA consisted of small RNAs (<200nt), with 
minor amounts of 18S and 28S rRNA (Figure lA). 
Next, small (<200nt) RNAs were isolated from the 
100000^ sedimented material, and processed for analysis 
by NGS. To assess whether selective RNAs were released 
into the extracellular space, smaU RNA from ceUs that 
had released the shuttle RNA was analyzed in parallel. 
Strand specific small RNA sequencing libraries were 
prepared and pre-amplified as described (36). Fragments 
with inserts between 15 and 70 nt were gel purified, after 
which the first 35 bases of the hbrary were read with a 
single read, with the lllumina platform. Fragments >70nt 
were not included to prevent the large number of tRNAs 
(~70nt in size) present in the cellular RNA fraction to 
dominate the sequencing reaction, which would hamper 
quantitative comparison of the cellular and shuttle RNA 
fractions. Sequencing data were processed using a 
modified version of the E-miR pipehne (36). The se- 
quences were mapped to the mouse genome (NCBI137/ 
mm9), using bowtie allowing for two mismatches and 
three alignments per sequence read, foUowed by a 
crossmapping correction (37). 

To vahdate the robustness of our sequencing results, we 
replicated the presented experiment with a slightly differ- 
ent sample preparation method (including an extra RNA 
fragmentation step). We found strong concordance when 
ranking the different RNAs based on their abundance in 
the shuttle RNA. We calculated the correlation between 
the (ranked) small shuttle RNA sequencing data of the 
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Figure 1. Shuttle RNA mainly consists of small RNA species and contains relatively low amounts of miRNA. Cell culture supernatant of DC-T cell 
co-cultures was subjected to several differential centrifugation steps, followed by RNA isolation from the lOOOOOg sedimenting material (shuttle 
RNA). (A) Bioanalyzer electropherogram of total shuttle RNA derived from DC-T cell co-cultures. The composition of small (<200nt) shuttle and 
cellular RNA was further analyzed by deep sequencing. (B) Sequence reads uniquely aligned to the genome were annotated to known non-coding 
RNA transcripts as annotated in Ensembl using the E-miR pipeline. Data represent the percentages of reads in the indicated categories of annotated 
transcripts in shuttle RNA (grey bars) versus cellular RNA (black bars) calculated relative to the total sum of reads for Ensembl non-coding RNA 
transcripts. (C) Distribution of shuttle and cellular RNA transcripts over the different RNA biotypes within the Ensembl-annotated miscellaneous 
RNA (misc_RNA) category. Data are represented as the percentages of annotated transcripts as in (B). 



replicates. A high correlation was observed between the 
replicates (Pearson's correlation 0.97, P-value <2.2e-16; 
Supplementary Figure SI). 

miRNAs are underrepresented in shuttle RNA compared 
with cellular RNA 

An overview of the sequencing data analysis and results 
can be found in Figure 2. We found that 8.3 x lOE* 
cellular and 27.5 x lOE'' shuttle RNA sequences aligned 
to the mouse genome (Figure 2). Uniquely aligned 
sequence reads were annotated to known non-coding 
small RNA transcripts as annotated in Ensembl. 
Remarkably, the percentage of sequence reads mapping 
to this category of RNAs was much lower in shuttle 
RNA released from cells (0.5%) compared with cellular 
RNA (31%). Within this category of Ensembl-annotated 
non-coding small RNAs, we analyzed the distribution of 



RNAs over the different biotypes. We found that the 
majority of sequences present in the cellular small RNA 
population represented miRNAs (Figure IB), which is in 
agreement with previous studies (36,41). In contrast, the 
percentage of miRNAs in shuttle RNA was very low. 
Instead, this fraction contained a relatively higher percent- 
age of small ribosomal RNA (rRNA) and RNA belonging 
to the miscellaneous category (misc_RNA) (Figure IB). 
Within this category of misc_RNA, we found that struc- 
tural SRP-RNA, vRNA and Y-RNA were most abundant 
in shuttle RNA (Figure IC). 

The miRNA distribution in shuttle and 
cellular RNA is different 

Since the initial discovery of miRNAs in ceU-derived 
vesicles (17), many research groups have studied the oc- 
currence of this type of regulatory RNAs in vesicles 
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Truncated sequences 

Cells: 9,887,963 
Shuttle: 28,820,290 



Aligned sequences 

Cells: 8,273,374 (83.7%) 
Shuttle: 27,500,560 (95.4%) 




Annotation exons, introns, and repeats 



EnsembI small 
non-coding RNA 

Cells: 2,568,765 (31 .0%) 
Shuttle: 134,932 (0.49%) 



Not EnsembI small 
non-coding RNA 

Cells: 5,704,609 (69.0%) 
Shuttle: 27,365,628 (99.5%) 



Aligned and in regions with > 5 
reads 

Cells: 4,880,148 (59.0%) 
Shuttle: 18,063,486 (65.7%) 



miRNA 

Cells: 1,471,703 (57.3%) 
Shuttle:1 0,784 (7.8%) 



not mlRNA 

Cells: 1,097,062 (42.7%) 
Shuttle:124,148 (92.2%) 



Exons 

Cells: 1,967,028 (40.3%) 
Shuttle: 1,255,824 (6.7%) 



Introns 

Cells: 512,776 (10.5%) 
Shuttle: 3,502,030 (19.4%) 



Repeats 

Cells: 1,453,876 (29.8%) 
Shuttle: 4,946,81 1 (27.4%) 



Exons of protein 
coding regions 

Cells: 967,229 (49.2%) 
Shuttle: 1,060,780 (84.5%) 



Exons of non-protein 
coding regions 

Cells: 999,799 (50.8%,) 
Shuttle: 195,044 (15.5%) 



Figure 2. Overview of the sequencing data analysis and results. Indicated are the number of reads in each analysis step for small cellular versus 
shuttle RNA and the relative amount of reads calculated as a percentage of the reads detected in the previous steps. 



derived from various cell types (as listed in the ExoCarta 
database (42)). Although miRNAs formed only a small 
minority within the total pool of shuttle RNA sequences 
(Figure IB), we examined whether a specific set of 
miRNAs was released during DC-T cell interactions. 
The miRNAs occurring in cells of the DC-T cell co-culture 
and shuttle RNA released by these cells were ranked based 
on the E-miR data and the 10 most prevalent miRNAs in 
each population were compared (Table 1). miRNAs were 
named according to the 5' (5p) or 3' (3p) arm of the 
hairpin where they originated from (36). Certain 
miRNAs, such as miR-29a and miR-31, were abundant 
in both cellular and shuttle RNA. However, several highly 
abundant cellular miRNAs, such as miR92a-l and let-7b, 
occurred only in very low abundance in shuttle RNA. 
Conversely, the highly abundant shuttle miRNAs 
miR-223, miR-142 and miR-93 were less abundant in 
cellular RNA. These data indicate that dissemination of 
miRNAs into the extracellular space is not a random 



process, but that cells actively sort selective miRNAs for 
extracellular destination. 

Shuttle and cellular RNA contain different sets of 
transcripts 

Since the shuttle RNA population contained only 0.5% 
Ensembl-annotated smaU non-coding RNA transcripts, 
we next aimed to define which other types of transcripts 
were present in this population. By comparison with 
genomic locations obtained from EnsembI Biomart (39), 
we found that 26.1% of the shuttle RNA reads (against 
49.9% in cellular small RNA) mapped to introns and 
exons (Figure 2). Cellular small RNAs showed a higher 
frequency of exonic localization (including a large number 
of exons coding for miRNAs). However, the relative 
amount of exonic sequences mapping to protein coding 
regions was higher in shuttle RNA. In addition, shuttle 
RNA contained relatively more intronic sequences 
(Figure 2). Data have recently accumulated demonstrating 
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Table 1. Enriched miRNAs in cellular and shuttle RNA 



miRNA 



miRNA arm 



Abundant cellular miRNAs 



Rank 
in cells 



Rank in 
shuttle RNA 



miR-155 
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29 


miR-29a 


3p 
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1 


miR-92a-l 


3p 
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miR-31 
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5 


miR-15b 
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46 


miR-744 
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47 


miR-let-7b 
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119 


miR-191 
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miR-24-2 
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10 


miR-24-1 
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10 
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Abundant shu 


ttle miRNAs 






miR-29a 


3p 


2 


1 


miR-21 


5p 


15 


2 


miR-223 


3p 


41 


3 


miR-191 


5p 


8 


4 


miR-31 


5p 


4 


5 


miR-142 


5p 


47 


6 


miR-93 


5p 


48 


7 


miR-I03-l 


3p 


18 


8 


miR-103-2 


3p 


14 


9 


miR-24-2 


3p 


9 


10 



miRNAs were ranked according to the frequency of their occurrence 
(number of reads). Indicated are the 10 most prevalent miRNAs in cells 
or shuttle RNA and the corresponding ranking in the other sample. 



that protein-coding genes can also produce a complex set 
of non-coding RNAs (43), including transcripts 
originating from introns and from the untranslated 
regions (UTRs) (30,43). Interestingly, we observed that 
the percentage of reads in protein coding loci that 
aligned to 3'- and 5'-UTRs was much higher in shuttle 
RNA compared with cellular RNA (Figure 3A). Most 
of these transcripts had the same directionahty as the 
coding mRNA. Although their exact function is 
unknown, these UTR-derived small RNAs have been sug- 
gested to play a regulatory role in the attenuation or regu- 
lation of translation (44). 

Further analysis of the exonic sequences in shuttle and 
cellular RNA revealed that the large majority (84%) of 
exonic reads in shuttle RNA annotated to protein 
coding transcripts (Figure 3B), preferentially locating in 
the UTR regions of those transcripts (Figure 3A). Also the 
relative abundance of the small non-coding vRNA, 
Y-RNA and SRP-RNA was higher in shuttle RNA 
compared with cellular RNA. In contrast, lincRNA and 
miRNA are less abundant in shuttle RNA compared with 
cellular RNA (Figure 3B). Taken together, these results 
suggest that shuttle RNA is enriched in non-coding RNAs 
other than miRNAs and hncRNAs and that many of these 
are stiU to be annotated. 
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Figure 3. Exonic hits in shuttle RNA preferentially locate to UTR regions. (A) Exonic hits were analyzed for overlap with UTR regions. Indicated is 
the percentage of reads in shuttle RNA (left) and cellular RNA (right) mapping to the sense or antisense strand of UTR-regions, calculated relative 
to the total number of reads in Ensembl annotated exons. (B) The identity of exonic hits in shuttle (black) and cellular (grey) RNA. Indicated are the 
percentages of reads in the different RNA biotypes calculated relative to the total number of exonic reads. 
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Figure 4. Different distribution of RNA repeat sequences in shuttle and cellular RNA. Abundance and classification of sequencing reads that match 
different types of RNA repeats. Indicated are the percentages of all reads (A) or reads in the top 1000 of most abundant regions (B) in the indicated 
repeat types, as detected in shuttle (black) and cellular (grey) RNA. 



The distribution of RNA repeat sequences is different in 
shuttle and cellular RNA 

Since the majority of small shuttle RNA transcripts did 
not correspond to Ensembl-annotated small non-coding 
or protein-coding RNAs, we next investigated whether 
the as yet undefined transcripts corresponded to 
repeat-derived RNAs. The relative abundance of repeat 
RNAs in shuttle and cellular RNA was approximately 
the same (+/— 25%, Figure 4A). However, the distribu- 
tion over the type of repeat sequences in these two RNA 
populations was different. Shuttle RNA contained rela- 
tively large numbers of Long Interspersed Elements 
(LINEs), Long Terminal Repeats (LTRs) and simple 
repeat sequences, whereas the majority of repeat se- 
quences in the cellular RNA consisted of tRNAs 
(Figure 4A). To further classify the most abundant 
repeat-associated RNA species, we next considered only 
the top 1000 of most prominently present shuttle and 
cellular RNA. The differences in the distribution of the 
most abundant shuttle and cellular RNA were much more 
pronounced. Almost 90% of the 1000 most abundant se- 
quences in shuttle RNA (versus 33% in cells) consisted of 



repeat sequences (Figure 4B). The most prevalent types of 
repeats in shuttle RNA were tRNAs, LTRs, LINEs and 
simple repeats. LTRs, LINEs and simple repeats were 
relatively much more abundant in shuttle RNA 
compared with cellular RNA (Figure 4B). 

Specific tRNA fragments are abundantly and selectively 
present in shuttle RNA 

To investigate in greater detail the most abundant se- 
quences in shuttle RNA, all unique regions were ranked 
by abundance (number of reads found per region) and the 
RNA biotypes of the top 75 highest ranked hits were 
analyzed. For comparison, a similar top 75 ranking was 
made for the sequences found in cellular RNA (Table 2). 
We observed that specific snoRNA and niiRNAs were 
among the most abundant sequences in the cellular 
small RNA fraction, but not in shuttle RNA. However, 
as expected from previous data (Figure 4B), the shuttle 
RNA was more abundant in different types of repeat se- 
quences. The high abundance of tRNA hits in both the 
shuttle and cellular RNA fractions was remarkable, since 
we restricted our sequencing analysis to <70nt RNAs. 
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Consequently, the sequence reads were not likely to rep- 
resent mature full-length tRNAs. We next explored 
whether these reads could represent tRNA fragments 
(tRFs), which have recently gained interest due to their 
suspected regulatory nature (45). tRNA fragmentation is 
a specific process, with the composition, abundance and 
cleavage site varying per organism and cell type. 
Fragments of 18-22 nt, 30-35 nt and ~50nt locating at 
the 3'- or 5'- end of mature tRNAs have previously been 
described (45^7). In cellular RNA, we observed frag- 
ments located at the 5' or 3'end of the mature tRNA 
(Table 3). However, the most abundant tRNA hits in 
shuttle RNA were all located at the 5'end of mature 
tRNAs. Differences were also observed in the type of 
tRNAs occurring in the top 75 most abundant hits 
of cellular and shuttle RNA; whereas fragments of 
tRNA-Lys-AAA were abundant in both cellular and 



Table 2. Most abundant biotypes in cellular and shuttle RNA 

RNA biotype No. of reads No. of regions 



Cellular RNA 



snoRNA 


615 728 


20 


miRNA 


599 709 


20 


tRNA repeat 


586 583 


27 


Unclassified 


379 199 


4 


rRNA 


107 855 


2 


Protein coding 


103 877 


1 


vRNA 


15 776 


1 


Shuttle RNA 






tRNA repeat 


670 645 


47 


Simple repeat 


484945 


6 


LINE 


325 549 


5 


rRNA 


120 797 


5 


vRNA 


45 546 


1 


Protein coding 


26492 


4 


SRP-RNA 


18553 


2 


Y-RNA 


5929 


1 



Regions were ranked by abundance (number of reads per region) and 
the top 75 highest ranked hits were categorized based on RNA biotype. 
Indicated are the collective number of reads mapping to the different 
biotypes and the number of different regions over which these reads 
were distributed. 



shuttle RNA, fragments of tRNA-Asp-GAY were only 
abundant in cellular RNA. We observed that the 
majority of abundant tRNA hits in shuttle RNA 
covered larger regions of about 40-50 nt (Table 3 and 
Figure 5), caused by the presence of reads representing 
tRNA fragments (rather than full-length tRNAs) 
missing the first 5-1 5 nt of the mature tRNA. 
Remarkably, at the same genomic locations, the 
coverage in cellular RNA was restricted to ~30-35 nt frag- 
ments (Figure 5), indicating the presence of two different 
fragments of the same tRNA, of which one is uniquely 
present in shuttle RNA. 

A specifically cleaved part of vRNA is abundantly present 
in shuttle RNA 

Apart from the repeat and tRNA fragments, three other 
RNA biotypes were abundantly present in shuttle RNA 
(Table 2): (i) vRNA; (ii) SRP-RNA (7SL-RNA); and (iii) 
Y-RNA. All three RNAs are non-coding polymerase 111 
RNA transcripts. Due to our <70nt size restriction, we 
did not expect to find fuU-length transcripts of these three 
RNA biotypes. vRNAs are a family of RNAs found 
associated with the vault ribonucleoprotein complex 
located in the cytoplasm (48). The function of these 
vault particles are largely unknown, but they are 
thought to play a role in transportation of molecules, 
such as niRNA, from the nucleus to the cytoplasm and 
in drug metabohsm (e.g. in tumor cells). Mouse vRNA is 
141 nt long, but it was recently discovered that smaU vault 
RNAs (svRNAs) can be generated from vault non-coding 
RNAs through a DICER-dependent and 
DROSHA-independent mechanism (32). These svRNAs 
can downregulate expression of specific genes, similar to 
miRNAs. Interestingly, we found that the predominant 
coverage on the vRNA in exosomes was in only one of 
the internal stem loop structures (Figure 6A, B). In 
contrast, the coverage in cellular RNA was predominantly 
at the 3' and 5' ends of vRN A, resembhng the localization 
of the described human svRNAs (32). These data indicate 
that a specifically cleaved part of the vRNA is preferen- 
tially shuttled into the extracellular space. 



Table 3. Most abundant tRNA fragments in cellular and shuttle RNA 



tRNA fragment Location No. of reads No. of regions Region length (nt) 



Abundant cellular tRNA hits 



tRNA-Asp-GAY 


3'-end 


389 792 


14 


35 


tRNA-Lys-AAG 


5'-end 


106551 


6 


30-35 


tRNA-Lys-AAA 


5'-end 


36 995 


3 


30 


tRNA-Lys-AAG 


3'-end 


18 079 


1 


35 nt 


Abundant shuttle tRNA hits 










tRNA-Lys-AAA 


5'-end 


246 384 


6 


35 


tRNA-Lys-AAG 


5'-end 


134 033 


9 


50 


tRNA-Gly-GGA 


5'-end 


85 930 


7 


50 nt 


tRNA-Gly-GGY 


5'-end 


81 737 


9 


40 


tRNA-Val-GTG 


5'-end 


44160 


5 


50 


tRNA-Glu-GAG 


5'-end 


27 900 


5 


40 


tRNA-Val-GTA 


5'-end 


13 198 


2 


40-50 



Hits in tRNAs in the top 75 highest ranked hits were annotated with anticodon, location of the region within the tRNA coding sequence (3'- or 
5'-end), collective number of reads (accumulative for different copies of the same type of tRNA), the number of regions over which these reads were 
distributed, and the length of the regions. 
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Figure 5. Different coverage of tRNA genes in shuttle and cellular RNA. (A) Screen shots were taken from the UCSC genome browser (http:// 
genome.ucsc.edu/) and show examples of tRNAs. Custom tracks were added to show sequence coverage in shuttle and cellular RNA. An example is 
shown of a tRNA gene (tRNA-Lys-AAG) on which differential coverage was detected in shuttle RNA (top) and cellular RNA (bottom). The r-axis 
represents the coverage at each genomic position. (B) tRNAscan-SE-predicted (http://lowelab.ucsc.edu/tRNAscan-SE) secondary structure of the 
mouse tRNA gene (tRNA-Lys-AAG) presented in (A). Indicated are the regions with coverage in shuttle and cellular RNA (dotted line), the 
additional coverage in shuttle RNA (black line), and the anticodon (circle). 



SRP- and Y-RNA are highly enriched in cell-derived 
vesicles 

SRP-RNA (bound to the signal recognition particle) and 
Y-RNA (bound to the Ro ribonucleoprotein complex 
(49)) both function in intracellular transport of proteins 



and RNAs (50). We found a 28 nt Y-RNA fragment and 
26-50 nt SRP-RNA covering regions in shuttle RNA. The 
abundance of SRP- and Y-RNA in material released by 
cells into the extracellular space is remarkable, since these 
RNAs belong to a small group of host RNAs that is also 
selectively incorporated together with viral genomes into 
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Figure 6. Different coverage of vRNA in shuttle and cellular RNA. (A) Screenshot from the UCSC genome browser showing the sequence coverage 
on the vRNA gene in shuttle RNA (top) and cellular RNA (bottom). The ;'-axis represents the coverage at each genomic position. (B) 
MFOLD-predicted (http://mfold.rna.albany.edu/) secondary structure of mouse vRNA at 21°C. Indicated are the regions with predominant 
coverage in shuttle RNA (black line) and cellular RNA (dotted lines). 



the capsids of several different viruses (51,52). To establish 
a link between viruses and cell-derived vesicles regarding 
incorporation of these RNAs, we investigated whether 
full-length SRP- and Y-RNA released in our DC-T cell 
co-cultures were indeed associated to membrane vesicles. 
Hereto, the lOOOOOg material that was sedimented from 
the culture supernatant was further fractionated by 
gradient density ultracentrifugation. This allowed separ- 
ation of RNA associated to large protein complexes, 
which are retained in the bottom of the gradient, and 
RNA enclosed in membrane vesicles, which float to low- 
buoyant densities (1.12-1.18 g/ml (35); Figure 7A). 
Subsequently, the presence of SRP- and Y-RNA in the 
different density fractions was analyzed by RT-qPCR in 
independent experiments. For comparison, the presence of 
these RNAs was analyzed in the cellular <200nt RNA 
fraction obtained from the DC-T cell co-culture from 
which the shuttle RNA was derived. Invariant endogenous 
controls (reference genes) for qPCR analysis of 
vesicle-enclosed RNA are unknown. As an alternative, 
all samples were normalized to the total input quantity 
of RNA. Various miRNAs, which were relatively 
abundant in cellular and/or shuttle RNA (Table 1), were 
analyzed in parallel as a control. We observed that full- 
length SRP- and Y-RNA indeed associated to low-density 
cell-derived vesicles and not to high-density complexes. 
Moreover, these RNAs were highly enriched in vesicles 
compared with cells (Figure 7B). In contrast, aU tested 
miRNAs were relatively more abundant in cells 
compared with cell-derived vesicles, confirming the 
sequencing data (Figure 7C). Conclusively, SRP- and 
Y-RNA are highly abundant in cell-derived vesicles, 
similar to the enrichment of these host RNAs in virions. 
These findings may point to an evolutionary conserved 
mechanism by which cellular RNA and viral genomes 
are selected and/or stabilized in membrane vesicles with 
extracellular destination or by which transferred RNA can 
function in cells targeted by these vesicles. 



DISCUSSION 

The complexity and diversity of the pools of extracellular 
RNA in cell-derived vesicles and vesicle-free complexes is 
daunting. Many research groups have focused on 
analyzing miRNAs in material released from a broad 
range of ceU types. However, the size distribution of 
shuttle RNA, as shown here and by others (53,54), 
extends well beyond the 20-23 nt size of miRNAs and is 
mostly in the range of ~20 to ~200 nt. This indicates that 
also other smaU RNA species are released by ceUs into the 
extracellular space. By deep sequencing of the <70nt 
fraction of this shuttle RNA, we found a large variety of 
small non-coding RNA species representing pervasive 
transcripts or cleavage products overlapping with 
protein coding regions, repeat sequences, or structural 
RNAs. We extended this analysis with detailed evaluation 
of the sequence coverage of several of the most abundant 
RNA biotypes in shuttle RNA. This led to the discovery 
that cellular and shuttle RNA can contain different small 
RNAs derived from the same non-coding RNA, as was 
the case for tRNA and vRNA. Although no conclusions 
can be drawn regarding the absolute concentration of dif- 
ferent RNA species, the data allow comparison of the 
relative amounts of RNA types in the pools of cellular 
and shuttle RNA. The unequal distribution of the 
detected RNA species over cellular and shuttle RNA, 
combined with increasing evidence for their role in gene 
regulation strongly suggest that ceUs specifically release 
these RNAs to modify the function of target cells. 

miRNA 

From both the sequencing and RT-qPCR data it became 
clear that miRNAs were underrepresented in shuttle RNA 
compared with cellular RNA (Figures IB and 7C). 
However, the sequencing data revealed that the miRNA 
composition in shuttle RNA is not a mere reflection of the 
cellular miRNA (Table 1), indicating that a specific set of 
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Figure 7. Culture supernatant of DC-T cell co-cultures was subjected to consecutive differential centrifugation steps. The lOOOOOg pelleted material 
was loaded at the bottom of a sucrose gradient, after which low-density vesicles were floated to equilibrium density by ultracentrifugation. 
Low-density fractions (1.12-1.18g/ml) and high-density (1.26-1.28 g/ml) bottom fractions of the gradient were collected (A). Small RNA 
(<200nt) was isolated from the lOOOOOg pelleted material present in these fractions. (B, C) RT-qPCR analysis of Y- and SRP-RNA (B) and 
different miRNAs (C) in small cellular RNA (cells). lOOOOOg pelleted vesicles floating to low-density fractions (low dens) and lOOOOOg pelleted 
non-floating high-density material (high dens). All samples were normalized to input quantity of RNA. Relative expression levels were calculated 
based on the Ct values. Indicated is the fold increase ± s.d. of PCR product in floating or non-floating material relative to cells (set to 1). 



miRNAs is selected for extracellular release. Evidence 
exists that miRNAs associated to cell-derived vesicles are 
fmictional upon transfer to target cells (17,18). Although 
the relative amounts of released miRNAs detected here 
may appear low, significant effects on target cell function- 
ing could be imposed due to efficient and specific cellular 
targeting of vesicles and the capacity of miRNAs to cause 
large-scale and long-term modulation of gene expression. 
Interestingly, 5 out of the 10 most abundant shuttle 
miRNAs have vahdated target genes that play important 
roles in immune regulation. miR-93, for example, targets 
Stat3, which is involved in regulation of T cell responses 
(55), and miR-223 targets Mef2c, which is necessary for 
the transcriptional activation of interleukins during per- 
ipheral T ceU activation (56). The sequencing and qPCR 
studies described here have been performed on total 
shuttle RNA obtained from DC-T cell co-culture super- 
natant. Future studies wiU reveal which miRNAs (and 
other small non-coding RNAs) are released by DC or by 
T cells and which target cells can be modified by this 
shuttle RNA. 



RNA repeats 

We found that repeat sequences were highly enriched in 
shuttle RNA compared with cells and that especially the 
LINE, LTR and simple repeats were more abundant in 
shuttle RNA. Due to the size restriction of RNAs 
selected for sequencing, the detected transcripts do not 
represent fuU-length LINE or LTR. Small RNAs derived 
from these repeats (repeat-associated smaU interference 
RNA) have previously been shown to originate from 
various scattered regions within repeats such as SINE, 
LINE and LTR (57,58). Interestingly, retrotransposon 
RNA was recently found enriched in microvesicles 
derived from tumor cells (59). It is currently unknown 
whether this retrotransposon RNA can be transferred to 
other cells and whether it can insert into the target 
cell genome. However, the dissemination and active tar- 
geting of retrotransposable elements or fragments 
thereof to other cells may be an effective strategy of 
cells to modify genes and regulate gene expression in 
other cells. 
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A large number of sequences detected in shuttle RNA 
mapped to tRNA loci. Full-length tRNAs were, however, 
excluded from analysis by size selection of RNAs smaller 
than 70 nt. Although we cannot exclude low-level contam- 
ination with full-length tRNA, most tRNA hits will rep- 
resent tRNA fragments (tRFs). Many recent publications 
have demonstrated that tRFs are not merely degradation 
products but are specific cleavage products with versatile 
functions, such as the inhibition of translation and 
guidance of other RNAs (60). We here demonstrated 
that the tRF compositions in cellular and shuttle RNA 
were different, indicating selectivity in released tRFs. 
Earlier evidence that tRFs can be released from cells 
and transported to distant cells comes from the field of 
plant physiology, where specific tRFs were identified in 
the phloem sap of several different plants (61). The 5' 
end located tRNA fragments observed in shuttle RNA 
could represent tRNA halves, which are produced by a 
single cleavage event in the anticodon loop that cuts the 
mature tRNA into one 5' and one 3' halve. It has been 
shown that these halves are often not present at equal 
quantities, differ in functions, such as the induction of 
stress granule assembly (62), and can be recruited to 
distinct cytoplasmic structures (63). Our observation that 
different cleavage products of the same RNAs distributed 
unequally over cellular and shuttle RNA further 
strengthens the concept that cells select specific RNAs 
for extracellular release. We also observed that the 
sequence coverage on particular tRNA genes was different 
in shuttle RNA compared with cellular RNA. This could 
indicate that specific tRNA cleavage processes are 
involved in the generation of tRFs for dissemination to 
the extracellular space. 

Infrastructural RNAs 

We observed a remarkable enrichment of sequences 
matching SRP-RNA and the less studied infrastructural 
RNAs vault- and Y-RNA in shuttle RNA. Similar to 
tRNAs, all of these structural RNAs are transcribed by 
Polymerase III and function in association with cytoplas- 
mic proteins. The full-length forms of these RNAs are 
longer than 70 nt, and only fragments were detected by 
sequencing. Small vRNA fragments (svRNA) covering 
the 5' and 3' end of the vRNA have been detected in 
human cells (32), and we detected similar fragments in 
cellular RNA derived from mouse DC-T cell co-cultures. 
Interestingly, in shuttle RNA we detected coverage on a 
different stem-loop structure. Future studies should reveal 
whether this specific fragment of vRNA can also be pro- 
cessed by Dicer and whether it can function in regulation 
of gene expression. Indications based on qPCR analysis 
suggest that not only fragments, but also full-length forms 
of SRP- and Y-RNA were present in shuttle RNA, more 
specifically, in low-density vesicles released by the cells. 
The relative high amounts of full-length SRP- and 
Y-RNA (Figure 6B) in cell-derived vesicles may explain 
the higher enrichment of SRP- and Y-RNA observed by 
qPCR in comparison with fragment sequencing. Y- and 
SRP-RNAs belong to a small group of 5-6 host RNAs 
that are also selectively incorporated together with viral 



genomes into capsids of several different viruses (51,52). It 
is suggested that these host RNAs are encapsidated to 
enhance virus assembly, virion stabihty and/or viral infect- 
ivity. Many enveloped viruses exploit the existing routes of 
membrane traffic to leave the host cell. A relationship 
between endosomal and/or plasma membrane routing of 
viruses and cell-derived vesicles can therefore be envi- 
sioned (64). We hypothesize that these RNAs could play 
a role in the specific sorting of regulatory RNAs into 
cell-derived vesicles. Alternatively, these RNAs could sta- 
bilize the RNA content of cell-derived vesicles or guide 
regulatory RNAs for efficient functioning upon release 
into the cytoplasm of target cells. 

The abovementioned RNA biotypes make up only 
around 50% of the total transcripts present in shuttle 
RNA. We found that many of the other transcripts 
located to intergenic regions. One example is a highly 
abundant transcript matching a region located 3kB 
upstream of the Tial gene, coding for T cell intracellular 
antigen 1, which is involved in apoptosis and mRNA 
sorting into stress granules. This non-annotated transcript 
locates to a highly conserved region and might be classi- 
fied as PROMPT (28). The function of PROMPs is largely 
unknown, but could involve positively or negatively 
influencing the expression of downstream located genes. 
Future studies on the large number of (non-)annotated 
intergenic elements enriched in shuttle RNA may 
uncover additional gene regulatory sequences that are 
transferred between cells. Other issues that need to be ad- 
dressed in the future include the deep sequencing of the 
larger, 70-200 nt fraction of small shuttle RNA and 
determining the cellular origin (DC or T cell) of the 
observed RNA species. The RNA in DC- and T 
cell-derived vesicles could for example be separately 
analyzed after absorption of vesicles released in DC-T 
cell co-cultures onto beads coated with cell-type specific 
antibodies. 

Since the relative amount of small cellular RNA se- 
quences that map to 18S and 28S rRNA is low, we 
expect that experimentally induced RNA degradation in 
our studies is limited. Although some of the abundant 
RNA species detected in shuttle RNA are cleavage 
products derived from mature non-coding RNAs, it is 
not likely that shuttle RNA represents non-selective 
disposal of RNAs that have been degraded inside cells. 
Non-selective disposal of degraded RNAs would lead to 
a much more similar distribution of the degradation 
products over cells and vesicles than indicated by our 
deep sequencing results. In fact, we even observed a 
highly unequal distribution of different degradation 
products of the same RNAs over cellular and shuttle 
RNA. Furthermore, evidence accumulates that (partly) 
degraded RNA fragments, such as tRFs, can also act as 
regulatory RNAs influencing gene expression. 
Uncontrolled release of these fragments would therefore 
impose a risk on disturbing cellular homeostasis. 

By classification and quantification of the deep 
sequencing data, we gained important information on 
the type of RNAs selected by immune cells for extraceUu- 
lar release. This study revealed that shuttle RNA contains 
a wealth of different non-coding small RNAs. Many of 
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the highly abundant small non-coding transcripts present 
in shuttle RNA are evolutionary well conserved and have 
previously been associated to gene regulatory functions. 
These findings allude to a wider range of biological effects 
that could be mediated by shuttle RNA than previously 
anticipated. Gene regulatory functions of shuttle RNA 
should therefore also be sought beyond the action of 
specific miRNAs and mRNAs. 
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